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Abstract —The operation of groups of heavy-duty vehicles 
(HDVs) at a small inter-vehicular distance (known as platoon) 
allows to lower the overall aerodynamic drag and, therefore, to 
reduce fuel consumption and greenhouse gas emissions. However, 
due to the large mass and limited engine power of HDVs, 
slopes have a significant impact on the feasible and optimal 
speed profiles that each vehicle can and should follow. Therefore 
maintaining a short inter-vehicular distance as required by 
platooning without coordination between vehicles can often result 
in inefficient or even unfeasible trajectories. In this paper we 
propose a two-layer control architecture for HDV platooning 
aimed to safely and fuel-efficiently coordinate the vehicles in 
the platoon. Here, the layers are responsible for the inclusion 
of preview information on road topography and the real-time 
control of the vehicles, respectively. Within this architecture, 
dynamic programming is used to compute the fuel-optimal speed 
profile for the entire platoon and a distributed model predictive 
control framework is developed for the real-time control of the 
vehicles. The effectiveness of the proposed controller is analyzed 
by means of simulations of several realistic scenarios that suggest 
a possible fuel saving of up to 12% for the follower vehicles 
compared to the use of standard platoon controllers. 


I. Introduction 


The transportation of goods has been fundamental to the 
world economic development and the demand for freight 
transportation, together with the global economy, is expected 
to increase in the coming years. However, the transport sector, 
due to the burning of fuel, is responsible for a significant 
amount of greenhouse gas and C0 2 emissions. In the European 
Union, the transport sector amounts to roughly 29% of the 
total C0 2 emissions and 11% of these emissions are directly 
accountable to road freight transportation |3j, [ 171. Globally, 
the C0 2 emissions linked to the surface (road and rail) freight 
transport sector are expected to increase up to 347% in the 


next 40 years if no measure is taken 1191. In order to contrast 
this increase and the related impact on the climate change, 
governments all over the world are agreeing in setting stringent 
limitations on greenhouse gas emissions connected to road 
freight transportation Q5). In order to cope with these 
limitations, heavy-duty vehicle (HDV) manufactures are facing 
numerous challenges. Furthermore, the expected increase of 
the oil price 00 and the need for maintaining competitiveness 
require them to design vehicles and technologies that are 
increasingly fuel-efficient. The fuel cost for an HDV fleet 
owner, in fact, accounts roughly for the 35% of the total 
cost of owning and operating a vehicle 0. Therefore even 
a reduction of a few percent of the fuel consumption would 
lead to significant saving. 
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An effective method to reduce fuel consumption and, con¬ 
sequently, greenhouse gas emissions, is HDV platooning. By 
operating groups of vehicles at small inter-vehicular distances, 
the overall aerodynamic drag can be reduced. As about one 
fourth of the HDV fuel consumption is related to aerodynamic 
drag |2()| , platooning can have a large effect on the fuel 
consumption. Indeed, experimental results in |5| and 1131 
have shown a reduction in fuel consumption up to 7%. 
However, in order to safely operate HDVs at the short inter- 
vehicular distances required for platooning, automation of the 
longitudinal dynamics is necessary. 

In this work we present a novel control architecture for 
fuel-efficient and safe HDV platooning. Vehicle platooning 
is not a new control problem. The first works on vehicle 
platooning appeared in the sixties, e.g., fl8) , |28| , p0| . The 
main focus of these early works was the theoretical study of 
the dynamics of a string of vehicles with local information, 
with a particular attention on the study of string stability, 
i.e., the attenuation of disturbances in position, speed and 
acceleration along the string of vehicles. The vehicle pla¬ 
tooning concept received the first application interest under 
the Partners for Advanced Transportation Technology (PATH) 
project (32| l, where platooning (of passenger vehicles) has been 
investigated as a means to increase highway throughput. Under 
this project a control architecture based on vehicle-to-vehicle 
communication has been developed for the platoon formation 
and maintenance of fully autonomous vehicles [22]|. Although 
the environmental aspect was not the focus of the project, 
noteworthy results on fuel reduction due to HDV platooning 
have been reported G3- Since then, new projects and related 
publications have appeared with focus on different aspects of 
HDV platooning such as congestion, safety, fuel-efficiency and 
user-acceptance 0lj, (33}. 

In the more recent COMPANION project |Jj, where this 
paper finds its place, the fuel-efficiency of HDV platooning 
is the main focus. The project is not limited to pursue the 
efficiency of a single platoon, but rather to create a complete 
fuel-efficient freight transportation system. This led to the 
development of a system architecture aimed to divide this 
complex problem into hierarchical solvable subproblems (5J, 
(27}. An adaptation of such architecture has three layers, 
namely the mission planner, the platoon controller and the 
low-level vehicle controller, defined as follows: the mission 
planner is responsible for the optimal routing of the HDVs and 
their synchronization in order to create and dissolve platoons 
in optimized meeting points. This problem has been addressed 
in |26| where the authors propose a distributed framework for 
the synchronization of single HDVs and platoons on the road 
network. The platoon controller of each platoon receives from 
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the mission planner the optimal route and the average speed 
per link that the platoon should track. Therefore it controls the 
vehicles’ dynamics and it computes the inputs for the low-level 
vehicle controllers of each HDV. In 0 a distributed control 
framework over the platoon suitable for HDV dynamics is 
presented. However the role of external factors, such as slopes, 
is not taken into account. 

Because of the large mass and the limited power of HDVs, 
altitude variations have a significant impact on their behavior. 
Even small slopes produce such large longitudinal forces on 
the HDVs that they are often not able to keep constant speed 
during uphill segments (because of limited engine power) and 
during downhill segments without applying brakes (because 
of the significant inertia). Hence it is common that HDVs 
have to brake and therefore waste energy in order not to 
overcome the speed limit during downhill sections. This has 
been addressed in pT) in order to design a control system that 
optimizes the fuel consumption of single HDVs driving over 
hilly roads. In this work the authors showed how, by using 
look-ahead control based on road topography information and 
speed limits, it is possible to reduce the fuel consumption 
of a single HDV up to 3.5%. Slopes become more critical 
in the case of HDVs driving in a platoon formation. In |8J, 
the authors point out how the existing look-ahead strategies 
for single HDVs are not necessarily suitable for a platoon 
and that a dedicated approach is required. This is due to 
the fact that the additional requirement of keeping a small 
inter-vehicular distance between vehicles collides with the 
fact that HDVs experience significantly different longitudinal 
forces (e.g., gravity force depending on their mass and current 
road slope and air drag resistance depending on the distance 
from the previous vehicle). This appeared evident in the 
experimental results of a three-vehicle platoon driving on a 
highway presented in 0 - Even though in this work the HDVs 
have similar characteristics, the author highlighted how the 
use of feedback controllers in particularly hilly sectors of the 
highway could lead to an increase of the fuel consumption 
of the follower vehicles compared to the case in which are 
driving alone. These experimental data are further analyzed in 
Section[II]in order to obtain a good understanding of the role of 
the road gradient on HDV platooning. This analysis provides 
a motivation for the development of a novel cooperative look¬ 
ahead control for HDV platooning with the specific objective 
of fuel-efficiency, for which some early results have been pub¬ 
lished in p4) . Hence, this leads to the following contributions 
of the current paper. 

First, a control architecture for the fuel-efficient and safe 
control of an HDV platoon is presented. The control architec¬ 
ture is divided into two layers, namely the platoon coordinator 
and the vehicle controller layers. The platoon coordinator 
computes the fuel-optimal speed profile for the entire platoon 
by taking into account information about the road ahead. This 
optimal profile is communicated to the decentralized vehicle 
controller layer that safely tracks it and computes the real-time 
inputs for each vehicle in the platoon. 

Second, two receding horizon strategies within this control 
architecture are developed. The platoon coordinator relies on a 
dynamic programming (DP) formulation J9] that exploits pre¬ 


view information on the road topography and speed limits to 
compute a speed trajectory defined over space that is safe and 
fuel-optimal for the whole platoon. Here we emphasize that 
the platoon coordinator can handle heterogeneous platoons in a 
systematic way. The vehicle controller layer, instead, is solved 
though a distributed model predictive control formulation 03 
that tracks the speed trajectory and a certain gap policy while 
guarantying fuel-efficiency and safety. More precisely, it is 
proved that with this architecture no collision will occur within 
the platoon when up to one vehicle is controlled manually. 

The performance of the proposed control architecture is 
finally evaluated through extensive simulations motivated by 
real experimental scenarios and comparisons with existing 
approaches for speed control and spacing policy are presented. 

The rest of the paper is organized as follows. In Section [II] 
we analyze the experimental results presented in (5). In 
Section III we present the vehicle and platoon models used 
in the controllers, whereas in Section [TV] we introduce the 
control architecture. The platoon coordinator and the vehicle 
controller layers are discussed in Sections [V] and [VI] while 
their performance is studied in Sections |VII| |VIIl| and |IX| 
by means of simulations. Finally, conclusions are stated in 
Section Ixl 


II. Motivating experiment 

In this section we analyze the experimental results presented 
in [5] in order to reach a good understanding of the impact of 
the road gradient on HDV platooning and motivate the need for 
the design of a look-ahead control framework for fuel-efficient 
HDV platooning. 

In this experiment a platoon of three similar HDVs (same 
powertrain and mass of 37.5, 38.4, 39.5 tons, respectively) is 
driven over a 45 km highway stretch between the Swedish 
cities of Mariefred and Eskilstuna. The topography for this 
road is displayed in Figure |T] where the red color highlights 
the uphill and downhill sections where the slope is too large 
for a nominal HDV (whose parameters are displayed in Table [T] 
in Section | VII| ) to maintain a constant speed of 22 m/s 
without braking or exceeding the engine power limit. For 
the considered road, the steep sections represent 23% of the 
total length. Overall, the follower vehicles, by platooning, 
manage to save 4.1% and 6.5%, respectively. However in [5| 
it is shown that the fuel-efficiency drops significantly in the 
road sector where the slope is more varying. In this study 
we therefore analyze the behavior of the first two vehicles 
while driving over the particularly hilly stretch highlighted 
in Figure |T] as Sector A for which an increase of the fuel 
consumption of the second HDV of 4% compared to the case 
of driving alone has been reported. The behavior is shown in 
Figure [2] The first vehicle tracks a reference speed of 21.5 
m/s using cruise control and it switches to braking mode only 
when the speed limit of 23.6 m/s is reached, while the second 
vehicle tracks a headway gap (a distance proportional to its 
speed) from the first vehicle and it switches to braking mode 
only when the headway gap reaches a certain threshold (refer 
to (5) for a complete characterization of the controllers). In the 
analyzed sector three critical segments highlighted in Figure [2] 
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are identified where the use of feedback controls shows its 
limitations. 

• Segment 1\ due to the steep downhill the first HDV is 
not able to maintain the reference speed and therefore 
it accelerates while coasting. The second vehicle, while 
trying to track the headway gap policy, follows the 
same behavior. However, due to the reduced experienced 
air resistance, during the downhill the second vehicle 
accelerates more than the first one and, when the critical 
headway gap is reached, it brakes. In this case the co¬ 
ordination between the accelerations of the two vehicles 
would have the potential to avoid the undesired braking. 

• Segment 2: the headway gap deviates significantly from 
the reference one, due to a large relative speed at the 
beginning of the uphill segment and a change of gear 
during the segment. The second vehicle, in order to 
reduce the headway gap error, significantly increases the 
relative speed. Once the critical headway gap is reached, 
it brakes strongly. In this case, the prediction of the 
vehicles behaviors would have allowed the second vehicle 
to reduce the relative speed before reaching the reference 
headway gap and, therefore, to avoid the undesired brak¬ 
ing. 

• Segment 3: here the second vehicle shows a more critical 
behavior compared to the first downhill. In fact, during 
downhills, the vehicles’ actuators work close to saturation 
(small fueling and small braking) which is not suitable for 
feedback controller. Therefore in Segment 3 the control 
state of the second vehicle continues to switch between 
fueling and braking modes. In this case, the use of 
a receding horizon framework would have allowed to 
predict correctly the vehicle behavior depending on the 
slope and, by taking into account the actuators’ saturation, 
therefore, to obtain a smoother behavior of the vehicle. 

The analysis of these experimental results provides a strong 
motivation for the development of a cooperative look-ahead 
control strategy for HDV platooning based on a receding 
horizon framework where the road gradient and the vehicles 
ahead can be explicitly taken into account. 

III. Modeling 

HDVs are complex systems with a large number of inter¬ 
acting dynamics. Due to their heavy load, the braking and 
powertrain systems of an HDV have to generate and transfer 
extremely high torques. In this section we first present the ve¬ 
hicle system architecture upon which the proposed controller is 
designed. Second, we introduce the model of the longitudinal 
dynamics of a single vehicle and a platoon with a particular 
focus on the components that play a significant role for the 
fuel consumption. Finally, we present the fuel model used to 
estimate the fuel consumption. 

A. Vehicle system architecture 

The functioning of an HDV is guaranteed by a large number 
of system units that communicate with each other through 
the controller area network (CAN) bus. A simplified control 
architecture of an HDV is shown in Figure [3] where only the 



Figure 2. Experiment results presented in |5] relative to the first two vehicles 
of a three vehicle platoon driving along the Sector A highlighted in Figure |Tj 
The first plot shows the road topography, whereas the second plot shows 
the speed of the two vehicles; the third plot shows the real and reference 
(according to a headway gap policy) between the vehicles; finally the forth and 
fifth plots shows respectively the normalized engine torque for both vehicles 
and the normalized braking force for the second vehicle (the braking action 
of the first vehicle is not available). For additional details, see j5). 



Figure 3. Simplified system vehicle architecture. 


system units that are of interest for our work are displayed 
|5J. A more detailed description of a complete vehicle system 
architecture is given in 1241. The communication with the 
outside world relies on the wireless sensor unit (WSU). This 
unit shares real-time information with the other vehicles within 
the platoon. The global positioning system (GPS) computes 
the absolute position of the vehicle, while the radar measures 
the distance and relative speed between the current vehicle 
and the preceding one. The real-time state information of 
the platoon coming from the WSU, the GPS and the radar 
are fused by the data processing unit and sent to the high- 
level controller. The high-level controller computes the desired 
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Figure 1. Road topography for the 45 km highway stretch between the Swedish cities of Mariefred and Eskilstuna. The red color highlights the uphill and 
downhill sections where the slope is too large for the considered HDVs to maintain a constant speed of 22 m/s without braking or exceeding the engine power 
limit. 


acceleration and a boolean variable defining which low level 
controller should track it. In particular, the brake management 
system (BMS) controls the braking actuators, while the engine 
and gear management systems (EMS and GMS) control the 
engine, the gearbox and the clutch to provide the requested 
acceleration. 


B. Vehicle and platoon model 

In this subsection we derive the model of the longitudinal 
dynamics of a single vehicle and the platoon that is then used 
in the controller formulation. Using Newton’s second law, the 
longitudinal dynamics of a single vehicle can be expressed by: 

Tril'Oj — Fe,i "b -^b ,i "b -fg ,i (ct(Si)) + F r i ~b b’d.i {Vi , dfj , ^ ^ 

Si = Vi, 


where Vi and s, are the states of vehicle i, respectively, the 
speed and the longitudinal position, m* is its mass and F c ( , 
Fb,i, Fg i, F I} i and fj, are the forces acting on the vehicle. 
We collect and Si in the state vector Xi = [ Vi Sj] T ). More 
specifically, F c , and F b , are the control inputs and repre¬ 
sent the forces generated by the powertrain and the braking 
system, respectively. The engine force F c , is characterized in 
Section III-C while the braking force F b , is assumed to be 
limited by the road friction and therefore bounded by 


- miVigp < F b4 < 0 , ( 2 ) 

where p and r/, denote the (positive) road friction coeffi¬ 
cient and the braking system efficiency, respectively. Next, 
F gj i(a(s,;)) is the force due to the gravity, modeled as 

F g:i {a(si)) = —m,:gsin(a(si)), (3) 

where g is the gravity acceleration and a(si) the road slope at 
position Si. The rolling resistance is represented by F r l and 
is modeled as 

Ft,i tt r iriig ? ( 4 ) 

where c r is the rolling coefficient. Finally Fi^(vi,di) is the 
aerodynamic drag, modeled as 

Fd }i {vi,di) = -\pA v C D (di)vt , ( 5 ) 

where p is the air density, A v is the cross-sectional area of 
the vehicle and Cd is the air drag coefficient. In order to take 
into account the influence of the inter-vehicular distance on the 
aerodynamic force that plays an essential role in platooning. 


s? 



Figure 4. Experimental data |23| | and regression curve of the normalized drag 
coefficient experienced by an HDV as function of the distance to the previous 
vehicle. 


the drag coefficient Cd is defined as a function of the distance 
to the previous vehicle dj. This dependence is modeled by 

C D (di) = C D ,o ( 1 - „ CD '" . ) , (6) 

V Cd, 2 + di J 

where the parameters Cd,i and C\>.i have been obtained by 
regressing the experimental data presented in J23) . The effect 
of the short inter-vehicular distance on the leading vehicle is 
neglected since it is smaller than one on the follower vehicles. 
The experimental data and the regression curve are displayed 
in Figure [4] 

Remark 1. In this work we have chosen to model the air drag 
coefficient on the basis of the experimental data presented in 
relative to the second vehicle of a two buses platoon 
driving at 80 km/h. In the literature reports on air drag 
coefficient or fuel consumption measures based on both real 
experiments ^,0 425]/ and fluid dynamics simulator 29 ^ 
are presented. They show a reduction of the air drag coefficient 
for short inter-vehicular distances. How the reduction relates 
to the inter-vehicular distance varies. This variability has been 
attributed to the weather condition (e.g, temperature, humidity 
or wind) and the shape of the vehicles. 

The model of a platoon of N v vehicles is given by the 
combination of equations {!]) for i = 1,..., N v and the distance 
definition: 

{ oo, if i = 1 , 

, (7) 

Si —i Si Li —ll ^ ^ £•) 

where l, denotes the length of vehicle i. 
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engine speed [rpm] 

Figure 5. BSFC map for a 400 hp engine regenerated from (HJ The plot 
shows the BSFC expressed in g/kWh as function of the engine speed and 
torque. The dotted lines represent equal power curves, while the blue thick 
line represents the collection of the fuel-optimal operation points for various 
generated powers. 

C. Fuel model 

The powertrain is a complex system composed by engine, 
clutch, gearbox and final gear that allows to transform the 
fuel’s energy into longitudinal force. In this subsection we 
derive a simple model of the powertrain that captures the in¬ 
trinsic relation between consumed fuel and generated traction 
force. In the model derivation we ignore transmission energy 
losses and the rotational inertia of the powertrain components 
because they are assumed to be negligible when compared to 
the vehicle mass. 

Engine performance is typically described by the brake spe¬ 
cific fuel consumption (BSFC), that defines the ratio between 
consumed fuel and generated energy for various operation 
points (i.e., engine speed and generated torque). In Figure [5] 
we show the BSFC map for an HDV engine of 400 hp |31| 7 
where the dotted lines represent the collection of operation 
points with equal generated power. This map can be easily 
converted in one that defines the fuel flow Si as function 
of the engine speed uii and the generated engine power Pi, 
i.e.. Si = <j>i(u>i, Pi). By assumption, the engine power Pj, 
passing through the clutch, the gearbox and the final gear is 
completely transferred to the wheels. The rotational speed, 
instead, changes between the transmission components and 
is finally transformed into longitudinal speed by the wheels. 
Ultimately, under the assumption of no longitudinal slip, the 
vehicle speed Vi can be defined as Vi = kigiUJi, where fcj is a 
constant gain and % is the gear ratio of the gearbox. Therefore 
the fuel flow can be expressed as a function of the speed v t , 
the traction force F c , and the gear ratio gi as 

St = 4>i( , U '' , Pc.iVt] ■ (8) 

\ki9i ) 

In order to be efficiently used in the control design, the fuel 



Figure 6. The plots show the optimal fuel flow and engine speed as function of 
the generated power. In the first plot we also display the fuel model expressed 
in m obtained by the regression of the raw data. 

model is further simplified by removing the dependence of 
the fuel flow S t on the gear ratio g, through the introduction 
of an additional assumption: the gear ratio can be changed 
continuously on a unlimited span and the gear management 
system chooses the most efficient gear ratio. Hence, we 
redefine the fuel model as 

Si = min <pi (uii, F e i Vi) = <^ op t(9) 

CJi 

The resulting curve <^ 0 pt,*(•) is depicted in Figure [6] and is 
linearly regressed in order to obtain the fuel model used in 
the controller design, defined by 

Si — Pl.iFe.iVi T P0,i- (10) 

From this analysis we can also obtain the bounds on the 
generated power that are independent from the engine speed 
and the gear ratio: 

Pmin,i — — Pmax,i■ (11) 

In Figure [6] the two fuel models in ([9]) and ( fTO} , and the 
correspondent optimal engine speed are displayed. We note 
that the approximation error is negligible. 

IV. Platoon controller architecture 
The system architecture for the look-ahead HDV platoon 
controller is shown in Figure [7] The mission planner suggests 
routes and platoon opportunities. The platoon coordinator 
supervises the platoon behavior exploiting road information. 
The vehicle controllers execute the speed profiles for the 
individual vehicles. 

The platoon coordinator layer exploits available information 
on the topography of the planned route to find a fuel-optimal 
speed profile for the entire platoon, while satisfying the 
average speed requirement provided by the mission planner. 
Hereby, in order to capture the dynamics of the road topog¬ 
raphy, it considers a horizon of several kilometers and takes 
the constraints of all vehicles in the platoon into account. As 
a result, it can be guaranteed that every vehicle in the platoon 
is able to track the required speed profile. A single speed 
trajectory is computed by the platoon coordinator, representing 
the speed of the platoon. However, when this speed profile is 
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Figure 7. System architecture for look-ahead HDV platooning. 

specified as a function of space (i.e., position on the road) 
and the inter-vehicle spacing is chosen according to a pure 
time delay, every individual vehicle in the platoon can track 
this single speed profile. It is remarked that this layer can 
typically operate in a receding horizon fashion, providing an 
updated speed profile roughly every 10 seconds or when the 
recalculation is needed due to a strong deviation from the 
original speed profile. Finally, as this layer is not safety-critical 
and not related to a specific vehicle, it can be implemented in 
any of the platooning vehicles or even in an off-board road 
unit. In Section [V] we present a DP approach to formulate and 
solve the stated problem. 

The vehicle controller is responsible for the real-time con¬ 
trol of each vehicle in the platoon and is aimed at tracking the 
desired speed profile as resulting from the platoon coordinator. 
It also exploits the communication between vehicles of the 
assumed trajectories to ensure the proper spacing strategy. 
This layer guarantees the safety of platooning operations to 
for instance avoid collisions between trucks. Because of the 
safety critical aspect, this layer is implemented in a distributed 
fashion in each vehicle of the platoon. More precisely, each 
vehicle controller runs in the block named high-level controller 
in the system vehicle architecture shown in Figure [3] In 
Section [VI] a distributed model predictive control approach for 
this problem is discussed. 

Figure [8] shows how the optimization problems in the pla¬ 
toon coordinator and the vehicle controllers interact, and their 
mathematical structure. Note how the platoon coordinator, in 
order to have a good prediction of the consumed fuel over 
the horizon, uses an accurate non-linear model of the vehicle, 
while the vehicle controller layer, in order to enable fast 
computation necessary for the real-time control of the vehicle, 
uses a linear vehicle model. 

V. Platoon coordinator 

The platoon coordinator is the higher layer of the platoon 
control architecture. It takes as inputs the average speed re¬ 



Figure 8. Optimal control problems solved in the platoon coordination and 
vehicle controllers. 


quirement v from the mission planner and the current vehicles 
state Xi(t) from their vehicle controllers. By exploiting the 
available information on the planned route (i.e., slope data 
a s and speed limits t’^ ax ), it generates a unique feasible and 
fuel-optimal speed profile v s ’*(-) defined over space for all 
the vehicles within the platoon (i.e., v\'*(z) = V s ’*(z) for 
i = j ..... N v , where 0 is the space variable). Furthermore, 
according to safety criteria, it specifies the time gaps r t , 
defined as the time delay between two consecutive vehicles 
passing through the same point, i.e.. 


Si{t) = Si- i(f - Tj). 


( 12 ) 


Note that this spacing policy is consistent with the requirement 
that all vehicles have to follow the same speed profile over 
space. This can be easily shown by computing the time 
derivative of the left hand side, 


dsj(t) 
dt 


= Vi (t) =v s i (s i (t )), 


and the right hand side of ( fl2] i, 

dSi-i(t - Ti) 


dt 


= Vi-i(t - Ti) 

= vt-^Si-lit - Ti)) = vl_ L (Si(t)), 


(13) 


(14) 


where u|(s) denote the speed of vehicle i at space s. In fact, 
by combining the time gap definition ( fl2| with ( fl3| l and ( fl4| , 
we obtain v s (s) = u-_ 1 (s). 

The coordinator layer is implemented using a DP framework 
0 that runs in closed-loop. The parameters that characterize 
the DP problem are the discretization space Asdp, the horizon 
length //dp and the refresh frequency /dp. We also define the 
horizon space length as Sqp = fT DP As DP . 
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In the coming subsections we introduce all the components 
of the DP formulation, i.e., the vehicle model, the constraints 
on the input and states and finally the cost function. 


A. Platoon model 


The platoon coordinator layer uses a discretized version of 
the vehicle model <(TJ, where the discretization is carried out 
in the space domain using the implicit Euler approximation. 
The discretized vehicle model is: 




l(Zk) - Vj(Zk-l) 
Asr 


—Fe,i( z k) + Fb,i(zk) 


(15a) 


l(zk) 


t s i(zk) - t s i(z k - 1 ) 

Asnp 


- mig[s\n{a{zk)) + c r \ 

- lpA v C D {d\(zk)){vl{zk))\ 

=1, d5b) 


where z k is the discretized space variable, v\{z k ), Fffzk), 
F^^Zk) and d\{z k ) are the speed, the engine and braking 
forces and the distance to the previous vehicle expressed as 
function of space, respectively. 

The advantage of using the space discretization is that, by 
relaxing the average speed requirement, there is no constraint 
depending on time. The relaxation is done by removing the 
average speed constraint and introducing instead travel time 
over the horizon in the cost function, hereby using an appro¬ 
priate weighting. This allows to ignore the time dynamics and 
therefore reduce significantly the computational complexity. 

A drawback of the space discretization is that the distance 
definition {7} cannot be expressed in the space domain. In¬ 
stead, the following approximated expression, as function of 
the current vehicle speed v\(z k ), is used: 


is introduced. The practical effect of this constraint is to reduce 
the search space of the dynamic programming algorithm to one 
dimension rather then the number of vehicles in the platoon, 
enabling fast computation. 


C. Cost function 

The objective of the platoon coordinator layer is to define 
the optimal speed profile that minimizes the fuel consumption 
of the whole platoon, while maintaining a certain average 
speed. This is done by defining the cost function as the 
weighted sum of two terms: a first term Jf(v s (zj),u\(zj)) for 
j = k ,..., k + T/dp — 1 and i = 1,..., N v representing the fuel 
amount consumed by the platoon and a second term J t (v s (zj)) 
for j = k .+ // [)P — 1 representing the travel time over 
the horizon, i.e.. 


Jdp(v s (^)),u s (^-)) = M v *( z j)> u i{zj))+PJt(v*{zj)), (20) 

where f represents a trade-off weighf] The term 
Jf(v s (zj),u s i (zj)) is computed by using the fuel model 
taking also into account a final term representing the 
kinematic energy of the platoon at the end of the horizon: 


M vS (zj),K(zj)) 


N v p—1 

E 

j=k 


=E E + ^ 


N v 


-E^m 


mi(v s (z h+ H DP - 1)) 2 


2=1 


d\{z k ) = v\{zk)Ti - Zj_i. 


(16) The term J t (v s (zj)) is obtained by using the time model (|15b|i: 


In the DP formulation we refer to ( fl5a| as 

vl(zk-i) = fF{vl( z k),ul(z k )), where u\(z k ) is the 
input vector defined as u\(z k ) = [F^(z k ), TjA(z fc )] T . 


Jt(v s (zj)) 


k-\-HQp— 1 

E 


j=k 


Asdp 

V s (Zj)' 


B. Model constraints 

The platoon model is constrained by introducing bounds on 
the input and the speed. 

1) Input constraints: According to Q and ([TT]i, the engine 
and braking forces are bounded by the following constraints: 

Praia,iVi(z k ) < F e ^(z k ) < Pmax,i v i ( z k ), {\1) 

-mipigp < Fffzk) < 0 . 

In the DP formulation, we refer to these constraints as 

K(zk ) £Ui (z k ). 

2) State constraints: In order to take into account the road 
speed limits, the speed is bounded by 

^min (z k ) < v\{z k ) < ^max (z k ). (18) 

We refer to this constraint as v\{z k ) G V s {z k ). 

Moreover, in order to require all the vehicle to follow the 
same speed profile, the constraint 

v i(z k ) = v s (z k ), i = l,...,N v . (19) 


D. Dynamic programming formulation 

We now have all the elements to formulate the DP problem 
solved in the platoon coordinator: 


min J DP (v s {zj),u s (zj)) 

(21a) 

U s (zj) 


subj. to vl(zj- 1 ) = flMiz^.uKzj)), 

(21b) 

u l(zj) e U-izj), 

(21c) 

Vi(Zj) = F(Zj) G V s (Zj), 

(2 Id) 

Zk = Si(f), 

(21e) 

V s (z k ) = Ui(t), 

(2 If) 


for j = k ,.... k + Hup—l, where the equations ( |21e[ ) and ( ]2 1 lj > 
represent the initial conditions of the DP formulation. 

1 Instead of the constraint on the average speed of Figure [h] the parameter 
ft is tuned to give the desired average time 









VI. Vehicle controller 

This section focuses on the distributed model predictive 
controllers running in the vehicle controller layer. 

Each vehicle controller runs locally. Vehicle i receives the 
optimal speed profile v s, *(-) and the time gap r, from the 
platoon coordinator and state information from the preceding 
vehicle. By tracking the optimal speed profile and gap policy 
requirement, and satisfying a safety constraint, it generates the 
optimal state and input trajectories, respectively x*(-\t) and 
a*(jt), and the desired acceleration a*it) for the vehicle low- 
level controllers. The parameters that characterize the MPC 
formulation are the discretization time A/mpc> the horizon 
steps number T/mpc- the refresh frequency /mpc and the length 
of the horizon defined as Tmpc = //mpc A/mpc- 

In the coming subsections we introduce all the components 
of the MPC formulation, i.e., the vehicle model, the constraints 
on the input and state, the safety constraint and finally the cost 
function. 


A. Vehicle model 

In the MPC formulation the vehicle is described by 

Xi(t j+ 1 \t k ) = Axi(tj\t k ) + Bai(tj\t k ), (22) 


where 


A = 


1 0 
AiMPC 1 


B = 


A£mpc 

0 


The variables Xi(tj\t k ) = \vi(tj\t k ) s i (f i |4)] T and oq(£j|4) 

denote the predicted state (speed and position) and control 
input (acceleration) trajectories of vehicle i associated to the 
update time t k , respectively. We also introduce three additional 
trajectories associated to each update time t k that will be used 
later in the MPC formulation: 

• the optimal state trajectory x*(tj\t k ), 

• the state reference trajectory Xi(tj\t k ), 

• the assumed state trajectory Xi(tj\t k ), 

for j = k....,k + //mpc — 1 and the corresponding in¬ 
put control trajectories defined likewise. While the pre¬ 
dicted and optimal trajectories are function of the opti¬ 
mization variable, the reference and assumed trajectories 
are pre-computed. More precisely the reference trajectories 
Xi(tj\t k ) = [ Vi(tj\t k ) Si(tj\t k )Y and di(tj\t k ) are computed 
from the reference trajectory v s ’*(-) and the current position 
s(t k ) of the vehicle. In particular, Si(tj\t k ) is defined recur¬ 
sively as 


Si (tj | t k ) 


Sz(tj), j = k, 

Si(tj-i\t k ) + A£ M pctZ s ’*(si(£j-i|4)), j > k , 


while Vi(tj\t k ) is defined as 


Vi{tj\t k ) = V S ’*(Si(tj\t k )). 


The control input reference trajectory a,i(tj\t k ) is defined as 
finite differences of Vi(tj\t k ), i.e., 

a-,(tj\ t k ) = i|4) Vi(tj |4))A/mpo 


The assumed state and control input trajectories are computed 
from the optimal and real trajectories of the vehicle as 


Xi(tj | tk ) 


Xi(tj), j < k, 

x*(tj\t k -x), k < j < k + H UPCl 


(23) 


and di(tj\t k ) likewise. As mentioned at the beginning of this 
section, each vehicle communicates the assumed trajectory 
Xi(tj\t k ) to the follower vehicle. In this case, the use of 
the optimal trajectory computed the previous step reflects the 
assumption of a maximum communication delay of At mpc- 


B. Input and model constraints 

In order to take into account the bounds on the braking 
force |2) and the engine power ( fTT| i, as done in the platoon 
coordinator layer, the control input a, is bounded by the 
following non-linear constraint: 




Fext(Xi, Si_i) 


<di< 


Pi. 


Pext(Xi , 


, (24) 


where F ext (xi,Si- 1 ) denotes the summation of the external 
forces acting on the vehicle and is defined as 


Fext(xi, Si- 1) = - TOig(sin(a(si)) + c r ) 

2 d(Si— 1 Si li)tJi . 


(25) 


The control input is additionally bounded by a soft constraint 
in order to allow braking only if necessary, i.e., when the 
safety constraint (see section VI-C| > is activated or the braking 
is required by the platoon coordinator. This is formulated as 
follows: 

ai + €i> min(a Cji , a»), e* > 0, (26) 


where e, is the softening variable and a CJ is the coasting 
acceleration (i.e., no braking and fuel injection) and is defined 

as: „ 

_ Pi,min f Fext(Xi, Si— i) (^,n\ 

ttc i — T • ‘) 

rriiVi TO* 


In the MPC formulation we refer to the constraint (24 »,< 25 
as CLi(tj\t, k ) £ Ai(xi(tj\t k )) and to the soft constraint (261,( 27]) 
aS 0>i (tj \ t k ) T €i (tj\t k ) £ Ae J i(Xi(tj\t k )). 

The speed is bounded according to the constraint ( fl8j ) a 

^min (tj |^fc) — Vmax (Si(tj\t k )). 


In the MPC formulation, we refer to this constraint as 

Vi(tj\t k ) e v(si(tj\t k )). 


C. Safety constraint 

The platoon is intended to operate on standard highways 
where other vehicles are present. The designed controller 
therefore should be able to cope with cases where the platoon 
behavior deviates from the predicted one because of internal 
disturbances (e.g., gear shift) or external events (e.g., related 
to the traffic situation or a vehicle cutting into the platoon). In 
this section we focus on the safety problem, leaving to further 
work the study of how such events should be handled (i.e., 
autonomously or switching to manual driving). 

The platoon is considered safe if, whatever a vehicle in the 
platoon does, there exists an input for all the follower vehicles 
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such that collision can be avoided. The safety of the platoon 
is guaranteed by ensuring that the state of each vehicle lies 
within a safety set and it is firstly studied by considering two 
adjacent vehicles and later extended to the entire platoon. In 
here we consider the following vehicle continuous dynamics: 


Xi — 


Vi 

Si 





(28) 


where Vi, 5, and a * are the speed, position and acceleration 
of vehicle i, respectively. 

Let us now focus on the dynamics of two adjacent vehicles 
described by 



Xi-l 

Xi 


F(xi-i,Xi,ai-i,ai) 


f(x i - 1 ,d i - 1 ) 


(29) 

where the acceleration of the current vehicle a,; is the 
control input, while the acceleration of the previous ve¬ 
hicle di -1 is the exogenous input that can be regarded 
as a disturbance. We also introduce the admissible set 
X = {Ci x]] T : Vi- 1 > 0, Vi > 0, Si -1 - s z > k-i} as the 
set of all admissible states, where l t denotes the length of 
vehicle i. In order to obtain a closed form of the safety set, the 
following conservative approximations of the the exogenous 
and control inputs are introduced: 


di -1 G .4 p (xi_i) 
di G A\xf) 


[0, u max j], 

[flmin,i) Umax,* 
[0j dma x,i]i 


if Vi- 1 > 0 , 
if Vi -1 = 0, 


(30a) 


if Vi > 0, 
if Vi = 0, 


(30b) 


where a min4 , a miM , a max i and a maXji are lower and upper 
bounds on the minimum and maximum possible accelerations 
of vehicle i, respectively. Such bounds are computed under 
reasonable assumptions on the vehicles and road properties, 
i.e., the vehicles’ speed is limited (0 < v, < u max ), the 
admissible vehicles’ weight is bounded (G M) and the road 
slope a is bounded (|a| < a max ). For example, the bounds 
a min j and a m j n .j can be computed as follows: 


—min,i 

= 

min 

®min. 


m , a, d), 

0 < 4 I<'U m ax 


^ 6 ^ max ,d >0 

^rnin,i 

= 

max 

^min. 

,i(v, 

m, a, d), 


0 ^'U^'Umax 


L: 6 ^ max , 0 


where 


amin,i(v,m,a,d) = -prjig - g sin(a) - c r - \pA v C D (d)v 2 

and —gr/ig represents the maximum braking capacity of 
vehicle i. Note that, due to the definition of the bounds and 
because of the dominance of the — gipg term in the definition 
of amin,i, the following inequalities hold: 

—min,i — L 0, (32a) 

—max,i — ttmax,z* (32b) 


In order to guarantee the safety of the subsystem ( |29| ), we 
should guarantee that the state [x]_ 1 xJ} T always lies in a 
safety set S included in X, for any admissible trajectory of 


Figure 9. Contour plot of the safety set boundary dS. The variable di denotes 
the distance between the two adjacent vehicles. 


the previous vehicle. We now define the safety set S C X, 
displayed in Figure [9j as 

S = {CiC : gj{xi-!,Xi) > 0, j = 1, ...,4}, (33) 


where 

„ f- - X - - , Cl vl 

gi \Xi— 1, Xi J — Si —i Si H—l ~ T ~; 

1 ^ a min ,i 

g'2(Xi-l , Xi) Si li— i, 

g 3 (Xi- 1 ,X i ) = Vi- 1 , 

94(Xi-l,Xi) = Vi 

(34) 

and we state the following result: 

Lemma 1. Given the dynamic system ( |29[ > and the constraints 
( |30a| l and ( |30b[ > on the exogenous and control inputs respec¬ 
tively, there exists a control law di = 4>([xJ_ 1 xf] T ) G Af{xi) 
such that for all [xJ_ 1 (to) xJ(to)] T G S and a,;_i G A p (xi-i), 
the condition [xf_ 1 (f) xf(f)] r G S holds for all t > to. In 
other words, S is a robust controlled invariant set Ell- 


Proof. By using Nagumo’s theorem for robust controlled 
invariant sets (121, the lemma can be proved by showing that 
for all [Cl-IT e dS (defined as the boundary of S) there 
exists an di G A f such that, for all «,_ | G A p , the relation 

VSj (*<_!, ®<) Fifti, Xi—\,di—\, df) f 0 (33) 


holds for all j such that gj(xi-±, xf) = 0. Because of the 
structure of the problem, the control input di is chosen as 


ttmin ,ii tf Vi > 0, 

0, if Vi = 0, 


(36) 


for any \xj_^ iJ] T G dS and a,_i G A v {xi-\). We organize 
the proof by considering the Ci C T G dS defined by the 
activation of each gj(xi-i,Xi) > 0: 

• for [x'l_ l xlY such that g\(xi-\,Xi) = 0, and 

gj{xi-i,Xi) > 0, for j = 2,3,4, 


(Xi- 

■1 j Q>i— 

- 1 5 ) 


di- 1 

j - 

1 — 

—min, 2 —1 . 

I Vi -1 - 

1 - 

di- 1 

C?i 

7' 

IV 


—min,z—l . 




^min,2 / 
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where the equality and inequality hold because of how 
di is defined by ( |36| ) and g 3 (xi--i , xf) > 0. 

• for [xJ^xJY such that g 2 (xi-\,Xi) = 0, and 

gj(xi-!,Xi) > 0, for j = 1,3,4, 

'Vg 2 {xi- 1 ,Xi) T F(xi- 1 ,Xi,di- 1 ,di) = v t -i - vg > 0 , 

where the inequality holds by noticing that the combi¬ 
nation of gi(xi-i,Xi) > 0, gi{xi-\,Xi) = 0 and the 
relation fl32a[ i gives v t -1 > (a miM /a min ,i)£i- 

• for \xl_xxlY such that g 3 {xi-x,xf) = 0, and 

gj(xi-i,Xi) > 0, for j = 1, 2,4, 

Vg 3 (xi- 1 ,Xi) T F(x i - ll Xi,di-i,di) = di-x > 0 , 

where the inequality holds because of ( |30a[ i. The same 
can be verified in a similar way for [xj_x a;J] T such that 
9 i{xi-x,Xi) = 0 and gj(xi-x,xf) > 0 for j = 1,2,3. 

□ 


The choice of the safety set guarantees that the follower 
vehicle can react to the emergency braking maneuver of 
its predecessor, such that both vehicles come to a standstill 
without colliding. We now extend the result in Lemma 1 to 
the safety of the whole platoon. More precisely, we proof 
that whatever a vehicle does, there exists an input for all the 
follower vehicles, such that collision can be avoided. This is 
formalized by the following theorem: 

Theorem 1. Consider a vehicle with index i 3 < N v and 
all its follower vehicles i £ I = {i 0 + l,...,iV,,} sat¬ 
isfying the dynamics in ( |28| ). Then, there exists a control 
law ai = (f>{Xi,Xi- 1 ) £ i £ I such that for all 

[iJ_i(fo) ^J(fo)] T £ S and di 0 £ A p (xi 0 ), the condition 
[ij_ 1 (f) iJ(i)] T £ <5 holds for all t > io and all i £ I. 

Proof The application of Lemma 1 for i = io + 1 proves 
the existence of an input di £ A { (xi) that ensures that 
[xj_i{t) iJXf)] 1 £ S for all t > t 3 . Then, by noting that 
A t (xi) C A p (xi) according to ( [32] ), it follows that di £ 
_4 p (.T,).The theorem is then proven by induction over the 
vehicle index, hereby repetitively applying Lemma 1. □ 


This result is adapted to the MPC formulation in order 
to guarantee the safety of the platoon. More precisely, each 
vehicle, knowing the assumed state trajectory of the vehicle 
ahead, can compute the safety set for its own predicted state. 
By taking into account that, according to the definition of the 
assumed state in ( |23j ), the real state of the previous vehicle is 
know with a one step delay, the safety set S translates to the 
following safety constraint on each follower vehicle state: 




V?{t j+ 1 \tk) 

2a„ 


— Si—l(tj—l \t k ) 


D f- lfa-lfa) 


(40) 


min 


li— 1} 


for i = 2, ...,N V . Note that for safety purpose only the safety 
constraint for j = k is necessary. In fact it guarantees that, 
if at the update time tj- the current state of each follower 
vehicle is safe, then it is going to be safe also at the update 
time L+i • However, the safety constraint for j > k gives 


optimal trajectories that are safe over the whole horizon and 
therefore produces a smoother and more fuel-efficient behavior 
of the platoon. In the MPC formulation, we refer to the safety 
constraint <|40]) as f sa ie{xi{tj+i\tk)) > 0. 


D. Cost function 

The objective of the vehicle controller layer is to follow the 
optimal trajectory and the gap policy requirement provided 
by the platoon coordinator layer. This can be formulated by 
introducing the following cost function: 


JY PC ( x i(£tk),a,i{-,tk),ei{-,tk)) 

k+H M pc — 1 

= ^ ^ || a; i(fj \tk) — Xi—i (tj—Ti |£fc)| |^q 

j—k 

+ \\Xi(tj\tk) — %i(tj\tk)\\(l_Q')Q 
-f | ( Xj (tj \ tk) Oi (tj 11 k ) 11 n 
+ ||ei(f/|ffc)||p j 


where 


0 = 


0, if * = 1, 

C> if i = 2, ...,N V 


(41) 


and Ti represents the discretized version of the time gap r, 
(i.e., Ti = (t i /AfjyipcJ). The parameters Q, R and £ £ [0, 1] 
can be chosen in order to have a good trade-off between ref¬ 
erence trajectory, gap policy tracking and actuators excitation. 
The weight P related to the softening-variable of the constraint 
( |26| ) is chosen relatively large such that only the activation 
of the safety constraint / S afe(*i(fj+i|ffe)) > 0 can require a 
significant braking force. 


E. Model predictive control formulation 

We now have all the elements to formulate the 
problem: 

min (x{(-, ffc), ffc), £»(') tk)) 

subj. to Xi(t j+1 \t k ) = Axi{tj\t k ) + Ba,i(tj\t k ), 
di (tj | t k ) £ A{ {X-i {t j\t k ) ), 

(li (tj \ t k ) 4 ‘ £i(fj\tk') £ Ae, i (Xi (t- j |f k ) ) , 
Vi(tj\t k ) £ V(si(tj\t k )), 
fse£e(,Xi (fj +1 |ffe)) > 0 , if * > 2 , 
t k ) > 0 , 

Xi{t k \tk) — Xi(tf 


MPC 

(42a) 

(42b) 

(42c) 

(42d) 

(42e) 

(42f) 

(42g) 

(42h) 


where j = k,k + 11 m'C — 1 and ( |42h[ > represents the initial 
condition of the MPC problem. For implementation purpose 
the state-dependent constraint set in ( |42c| ), ( |42d[ > and ( |42c[ ) 
will be replaced respectively by Ai{xi{tj\t k )), A e ^(xi(tj\t k )) 
and V(si(tj\t k ))- Taking into account that the safety constraint 
{42]]) is quadratic and convex, the MPC problem can be 
recasted into a quadratic constraint quadratic programming 
(QCQP) problem for which efficient solvers exist. 

The output of the vehicle controller is the desired acceler¬ 
ation a*(t k ) (defined as a*{t k ) = a*(t k \tk), where a*f\t k ) 
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is the optimal input trajectory resulting from the MPC) and a 
boolean variable ab r ,i defined as 

_ fl> if a i {tk) < /■«! 

bl,i \o, if a*(t k ) > a* ti (t k \t k ), 

that states if the desired acceleration should be tracked by the 
BMS or the EMS. 

VII. Performance analysis of the platoon 

COORDINATOR 

In this section we analyze the performance of the platoon 
coordinator (as presented in Section [V] and shown in Figures 
[7] and [8j by focusing on fuel-efficiency. We compare its 
performance with other standard controller setups. To make 
the analysis independent from the low-level tracking strategy, 
we assume in this section that the HDVs can follow exactly 
the speed trajectories and spacing policies defined by the high- 
level controllers. 

A. Experiment setup 

The comparison is done by using as benchmark the scenario 
introduced in Section [II] We therefore consider a platoon of 
two HDVs driving over the 45 km road stretch shown in 
Figure [T| and investigate the controller performance for both 
homogeneous and heterogeneous platoons. The performance 
metrics chosen to compare the different control configurations 
are the energy and the fuel consumed by the HDVs. In 
some comparisons the consumed energy is preferred over the 
consumed fuel because it can be directly related to the energies 
dissipated by the various forces (i.e., gravity, rolling, drag and 
braking forces). 

The control configurations considered in the comparisons 
include three control strategies and three gap policies. In detail, 
the following control strategies are considered: 

• cruise control (CC): the first vehicle keeps the constant 
reference speed vqc on low-grade slopes. If the uphill 
slope is too large to maintain constant speed, the en¬ 
gine generate the maximum power P max until the speed 
reaches Vqc again. If the downhill slope is too large 
to maintain constant speed without braking, the engine 
coasts (i.e., does not inject any fuel, generating therefore 
the minimum power P mm ) until the speed reaches vqc 
again. However, if the HDV reaches the speed limit t; max , 
the brakes are activated in order not to overcome it; 

• look-ahead control (LAC): the first vehicle exploits the 
slope information of the road ahead in order to minimize 
its own fuel consumption. 

• cooperative look-ahead control (CLAC): the first vehicle 
follows the speed profile generated by the platoon coor¬ 
dinator proposed in this paper. 

The following gap policies are considered: 

• space gap (SG): the second vehicle keeps a constant 
distance dsc. from the first vehicle; 

• headway gap (HG): the second vehicle keeps a constant 
headway time thg from the first vehicle, i.e., it keeps a 
distance proportional to its speed (dnc(t) = thg'gW); 


• time gap (TG): the second vehicle keeps a constant time 
gap ttg from the first vehicle according to ( | 1 2\ . 

In order to be able to maintain exactly the desired gap 
policies as previously assumed, the second vehicle is allowed 
to overcome the theoretical maximum engine power P max y 
and to brake if necessary. In addition, in order to obtain a fair 
comparison it is ensured, by tuning the trade-off parameter 
/? of the LAC and CLAC formulations (see ©)> that the 
different control strategies have the same average speed v and 
the parameters cisG- thg and t T g are chosen such that the vehi¬ 
cles in the different gap policies have the same distance when 
driving at constant speed v (i.e., d sg = vthg = vtjq — If). 
Finally in order to remove the influence of the residual 
kinematic energy, the initial and final speeds are constrained 
to be the same in all the controller configurations. 


B. Fuel-efficiency analysis for different control strategies 

In this section we present the results of the platoon behavior 
for the three different control strategies, while keeping a TG 
policy (ttg = 1.4 s). In the first part, as in the motivational 
example of Section [II] we focus on the homogeneous platoon 
scenario, while in the second part we consider two heteroge¬ 
neous platoons (i.e., where the second vehicle is respectively 
heavier and lighter than the leading one). 


Table I 

Vehicle’s paramters 


Parameter 

Value 

mass {rrii) 

40 t 

length (li) 

18 m 

roll coefficient (c r ) 

3 x 10 -3 

vehicle cross-sectional area ( A v ) 

10 m 2 

maximum engine power (P m ax,i) 

298 kW 

minimum engine power (Pmin i) 

-9 kW 


We now consider a platoon of two identical vehicles, 
whose parameters values are displayed in Table [I] We start 
the comparison by analyzing the comprehensive bar diagram 
displayed in Figure 10 representing the energy consumed by 
each vehicle of the platoon for the three control strategies (the 
corresponding fuel consumption is displayed in the central 
column of Table E- This energy is normalized respect to 
the energy consumed by a single vehicle driving alone using 
CC. The consumed energy is additionally split into various 
components representing the energy dissipated by each force, 
namely the gravity, roll, drag and braking force. We can first 
notice how the second vehicle, for all the control strategies, 
consumes less energy compared to the first one, due to the 
significant reduction of the energy associated to the drag force. 
Second, comparing the three control strategies, we can observe 
how the use of the LAC allows both vehicles to save energy, 
respectively 3.5% and 6.4% compared to the use of the CC. 
Instead, by switching from the LAC to the CLAC, the first 
vehicle consumes 0.1% more energy, while the second one 
saves 3.7% of energy; therefore the platoon, given by the 
union of the two vehicles, saves 3.6% of energy. This result 
is in line with our expectation since the LAC optimizes the 
fuel consumption of the first vehicle, while the CLAC targets 
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Figure 10. Comparison of the energy consumed by each vehicle of a platoon 
(mi = m 2 = 40 t), for the three control strategies, namely CC, LAC and 
CLAC, while keeping a TG policy, driving along the 45 km road displayed 
in Figure [T] Each bar represents the consumed energy normalized respect to 
the the energy consumed by a single vehicle driving alone using CC. The 
consumed energy is split into various components representing the energy 
dissipated by each force, namely the gravity (Eg), roll (E r ), drag ( K ( \) and 
braking (Eb) force. 


the reduction of the fuel consumption of the entire platoon. 
Consequently, the saving of the CLAC strategy with respect 
to the LAC strategy are expected to increase for platoons of 
more vehicles. Going into the details of the various consumed 
energy components, first we notice that the gravity and roll 
energy components are the same for both vehicles for all 
the considered control strategies. This is due to the fact that 
the gravity energy depends only on the difference of altitude 
between the initial and final points, while the roll energy only 
depends on the driven distance that is the same by experiment 
design specification. The drag energy, instead, is significantly 
different for the two vehicles because of its dependence on 
the distance to the preceding vehicle, while it is approximately 
the same for the different control strategies. What significantly 
changes between the different control strategies is the energy 
dissipated by braking. 

In order to understand the role of the control strategies in 
the braking usage in Figure [IT] we show part of the simulation 
results corresponding to the road highlighted as segment B in 
Figure |T] In this study we have chosen to focus on a downhill 
section because this is where the braking action is taking place. 
The comparison of the platoon behaviors follows: 

• CC: during the downhill, starting from speed vqc, the 
first HDV accelerates while coasting due to the large road 
grade. In the meantime the second vehicle has to brake 
slightly in order to maintain the time gap and compensate 
the reduced drag force compared to the first vehicle. At 
38.1 km, in order not to overcome the speed limit, both 
vehicles need to brake significantly; 

• LAC: by exploiting the topography information of the 
road ahead, the first vehicle reduces its speed before the 
downhill by anticipating the coasting phase such that the 
speed limit is reached only when the slope grade is small 
enough to stop accelerating while coasting and therefore 
it avoids braking. The second vehicle, as in the CC case. 



Figure 11. Comparison of the behavior of an homogeneous platoon (i.e., 
mi = m 2 = 40 t) for three different control strategies, namely CC, LAC 
and CLAC, while keeping a TG policy, driving along the Sector B displayed 
in Figure The first plot shows the road altitude, where the red color is 
used to highlight the sections too steep to keep a constant speed of 22 
m/s while respecting the power limit and avoiding braking; the second plot 
shows the speed profiles for the three control strategies followed by both 
vehicles (because of the SG policy); finally the third and forth plots show 
the summation between the generated power by the engine and the braking 
systems for the two vehicles and three control strategies; the black lines 
in such plots define the theoretical minimum and maximum engine power, 
respectively P m in,i and P max ,i (hence if the power crosses the lower power 
limit P m i n i, the respective vehicle is braking). 


has to brake slightly while the first vehicle is coasting but 
it also avoids the significant braking phase at the end of 
the downhill; 

• CLAC: since in this case the optimization is done consid¬ 
ering the fuel consumption of both vehicles, with respect 
to the LAC case the first vehicle starts to loose speed 
earlier before the downhill. This allows it to fuel sightly 
during the downhill, allowing the second vehicle to coast 
meanwhile and, as in the LAC case, to reach the speed 
limit only when the slope grade is small enough to stop 
accelerating while coasting. In this case both HDVs do 
not need to brake. 

Note that, in the case of longer downhill segments, the lower 
speed bound does not allow the vehicle to decrease the speed 
enough before the downhill in order not to hit the upper speed 
limit during the downhill. This is why in some sections of the 
45 km benchmark road, in the LAC case, the first vehicle and, 
in the CLAC case, both vehicles still need to brake. 

So far we have considered the case of an homogeneous 
platoon. What we want to investigate now is the the role of 
the different control strategies in the case of heterogeneous 
platoons. To answer this question, in Table [D] we have re¬ 
ported the normalized fuel consumption for the cases of two 
heterogeneous platoons and the same homogeneous platoon 
previously considered. More in detail, the HDVs have the 
same powertrain, but their masses vary between 35, 40 and 
45 t. Analyzing the table we can notice how in the case of a 
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heavier second vehicle the CLAC allows to save 10.8% of fuel 
compared to the CC, while, in the case of an lighter second 
vehicle, it allows to save 5.4%. However if we only analyze the 
last row we can note how, with the use of the CLAC, the order 
of the vehicles does not significantly change the normalized 
fuel consumption. 

Concluding, the proposed controller (CLAC) has a signifi¬ 
cant impact on the reduction of the energy and fuel consump¬ 
tion. In detail, the majority of the fuel saving is related to the 
reduction of energy dissipated by braking during the downhill 
sections. The impact of such a controller grows in the case of 
heavier follower vehicle. 


Table II 

Normalized fuel consumption of the vehicles in the platoon for different 
control strategies and scenarios (vehicle weights). The fuel is normalized 
respect to the fuel consumed by the respective HDV driving alone using 
CC. For the acronyms explanation refers to Section [Vll-A| [%]. 


mass 





a 2°d 

Pf lSt 

2 nd 

fl i st 

pi 2 nd 

35 t 

45 t 

40 t 

40 t 

45 t 

35 t 

CC 

100.0 

90.2 

100.0 

86.3 

100.0 

82.1 

LAC 

97.6 

84.9 

96.9 

80.6 

96.3 

77.2 

CLAC 

97.8 

78.0 

97.0 

77.4 

96.4 

76.7 


C. Fuel-efficiency analysis for different gap policies 

In the previous analysis we have always considered a TG 
policy. The aim of this section is to compare the platoon 
performance for different gap policies, namely the SG, HG and 
TG policies, while keeping the same control strategy (in the 
analysis we have considered CC). Note that in order to be able 
to follow the required gap policy the second vehicle is allowed 
to exceed the maximum engine power. In this section we only 
focus on the homogeneous platoon, since the results for an 
heterogeneous platoon are qualitatively the same. In Figure [I2| 
we show the comprehensive bar diagram representing the 
normalized energy consumed by each vehicle of the platoon 
for the three gap policies, while using CC as control strategy. 
Since the first vehicle uses the same control strategy, the 
energy consumption defers only for the second vehicle. It is 
interesting to notice that, similarly to the comparisons done in 
the previous section, the main difference between the energy 
consumption of the second vehicles is related to the energy 
dissipated by braking. More in detail the HG policy allows 
the second vehicles to save 1% over the SG policy, while the 
TG policy allows to save an additional 1.6% of energy. In 
order to understand the role of the gap policy on the braking 
energy, we show the platoon behavior driving over a synthetic 
hill composed by an uphill section with constant slope grade, 
a flat section and a downhill section with constant slope grade. 
The platoon behavior for such a hill is shown in Figure [13] 
Analyzing the second vehicle behavior for each gap policy, 
the following can be observed: 

• TG policy: as argued in Section [V] the time gap allows 
the vehicles to follow the same speed profile over space. 
That means that the generated forces and therefore the 


generated powers (because of the equal speed result) 
are equivalent except for a reduction of the air drag 
component in the second vehicle. Therefore the power 
generated by the second vehicle, as can be observed in 
Figure [13] is approximately a biased equivalent of that 
one generated by the first vehicle. 

SG policy: the space gap, instead, requires the vehicles 
to follow the same speed profile over time. An interest¬ 
ing consequence can be observed, for example, at the 
beginning of the uphill section shown in Figure 13 as 


soon as the first vehicle enters the uphill section and 
decelerates because of limited engine power, the second 
vehicle, which is still in the flat section, has to brake in 
order to respect the space gap requirement. In general, 
excluding the offset given by the drag power, every time 
the slope increases (in Figure [13] entering the uphill and 
leaving the downhill sections), the second vehicle has to 
generate less power than the first vehicle, while every 
time the slope decreases (in Figure [13] leaving the uphill 
and entering the downhill sections) the second vehicle 
has to generate more power than the first vehicle. As 
a consequence, the second vehicle has respectively to 
brake and to exceed the power limit in order to follow 
the required SG policy. 

• HG policy: the headway gap can be considered as a trade¬ 
off between a time gap and a space gap. In fact, for 
example, as soon as the first vehicle enters the uphill 
section and starts to decelerate, the distance between the 
two vehicles is allowed to decrease, but this decrease is 
not as fast as in the case of the time gap. 

The results obtained by the analysis of the platoon behavior 
in the case of the synthetic hill are valid also in the case of 
the original scenario. In conclusion, the time gap allows to 
save more energy compared to the space and headway gaps. 
In addition the time gap allows all the vehicles to follow the 
same speed trajectory in space and therefore it scales well with 
the number of vehicles in the platoon. The complete results 
for the normalized fuel consumption are reported in Table III 


Table III 

Normalized fuel consumption of the vehicles in the platoon for different 
control strategies and gap policies. The fuel is normalized respect to the fuel 
consumed by the respective HDV driving alone using CC. For the acronyms 
explanation refers to Section [VII-A| [%]. 



SG 

HG 


TG 










CC 

100.0 88.6 

100.0 87.7 

100.0 

86.3 

LAC 

96.9 82.7 

96.9 81.9 

96.9 

80.6 

CLAC 

97.0 80.4 

97.0 79.3 

97.0 

77.4 


VIII. Performance analysis of the vehicle 

CONTROLLER 

In this section we analyze the performance of the vehicle 
controller layer (as presented in Section [VI] and shown in 
Figures [7] and [8]) by focusing on the safety aspect. The analysis 
is based on the simulation result displayed in Figure 14 and 
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Figure 12. Comparison of the energy consumed by each vehicle of an 
homogeneous platoon (i.e., mi = m 2 = 40 t) for three different gap policies, 
namely space (SG), headway (HG) and time (TG) gap policies, while using 
CC as control strategy, driving along the 45 km road displayed in Figure [7] 
For the plots explanation refer to the caption of Figure |10| 
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Figure 13. Comparison of the behavior of a homogeneous platoon (i.e., 
mi = m 2 = 40 t) for three different gap policies, namely space (SG), 
headway (HG) and time (TG) gap policies, while using CC as control strategy, 
driving over a synthetic hill. For the plots explanation refer to the caption of 
Figure m note that the second plot shows only the speed trajectories of the 
second vehicle (the speed trajectory of the first vehicle coincides with that 
one of the second vehicle in the case of TG policy). 


Figure [15] where the leading HDV of a three vehicles platoon 
driving on a flat road brakes repeatedly with different braking 
profiles. Here we assume that the leading vehicle in the braking 
phases is manually driven and, therefore, the control system 
does not know a priori the braking profile. The considered 
vehicles are identical with the parameters as defined in Table [I] 

A. Safety analysis 

Here we focus on the safety analysis of the distributed 
vehicle controller layer and, in particular, we analyze the role 
of the safety constraint in various situations. In Figure [14] the 


Figure 14. Behavior of a three identical vehicles platoon driving on a flat 
road. The leading HDV brakes three times at 5, 25 and 55 s, with a braking 
deceleration of respectively 1, 2 and 3 m/s 2 for 0.9 s. The first plot shows 
the speed of the three vehicles: the second plot shows the distance between 
the vehicles and the respective safety distance computed using an adaptation 
of inequality (40) ; the third plot shows the summation between the generated 
power by the engine and the braking systems of the vehicles. 


leading vehicle is braking with deceleration of 1, 2 and 3 m/s 2 
for 0.9 s at respectively 5, 25 and 55 s. In the second plot of 
this figure, the effective distances and that ones that would 
activate the safety constraint (we will refer to it as the safety 
distance) are shown. First we can notice how, in line with our 
expectation, the second and third vehicles are braking (see 
the third plot) only when the effective distance touches the 
safety distance. In fact here we recall that, according to how 
the vehicle controller is designed (see Section VI-E| i, only the 
activation of the safety constraint or a braking request from the 
platoon coordinator can lead to a significant braking action. 
Consequently, during the first braking of 1 m/s 2 , both follower 
vehicles do not brake, despite the deviation of their states 
from the reference trajectories. During the second braking of 
2 m/s 2 , instead, the safety constraint of the second vehicle are 
activated and therefore it requires a braking action. Finally, 
during the third braking of 3 m/s 2 , the safety constraints of 
both follower vehicles activate and therefore they both brake. 
Note that the safety constraint is designed such that fuel- 
efficiency has priority on driver comfort. In fact, in this case, 
in order to be fuel-efficient, the braking action is required only 
when the platoon is in a safety critical situation. However, a 
priori knowledge of the braking profile of the first vehicle (e.g. 
having a model of the driver or handling the braking action 
autonomously) would have allowed to have a smoother and 
less intense braking action. 


In Figure 15 we consider a more challenging scenario in 
which the first vehicle brakes with higher intensity, simulating 
an emergency situation. More precisely it brakes at 5 s with 
a deceleration of 7 m/s 2 for 1 s and at 30 s with the same 
deceleration until it arrives to full-stop. We can notice how, 
also in this scenario the safety constraint in each vehicle 
controller layer activates the braking action and guarantees 
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Figure 15. Behavior of a three identical vehicles platoon driving on a flat 
road. The leading FiDV brakes a first time at 5 s for 1 s with a deceleration of 
7 m/s 2 and a second time at 25 s with a deceleration of 7 m/s 2 until arriving 
to full-stop. For the plot explanation refer to the caption of Figure |14| 




25 26 27 28 29 30 31 

space [km] 


no collision between the vehicles. 

IX. Performance analysis of the integrated 

SYSTEM 

In this section we analyze the simulation results displayed 
in Figure [16] of the platoon under the control of the integrated 
control architecture (i.e., platoon coordinator and vehicle con¬ 
troller). More precisely in this analysis we consider a platoon 
of three identical vehicles (whose parameters are defined in 
Table [I] driving along the Sector A highlighted in Figure [T] 
This is the same sector for which the experimental results in 
|5) are displayed in Figure [2] and analyzed in Section [II] 

At first glance, as expected from the platoon coordinator 
formulation, we can notice how all the vehicles follow the 
same speed and distance profiles in the space domain. Addi¬ 
tionally, in order to follow such profiles, we can observe in the 
last plot how the second and third vehicle, thanks to the air 
drag reduction, need to generate less power than the leading 
vehicle. We now continue the analysis by focusing on the three 
segments highlighted in Figure [16] 

• Segment 1: due to the steep downhill, all vehicles are not 
able to maintain the constant speed without braking and, 
therefore, accelerate. However the platoon coordinator 
requires the leading vehicle to fuel slightly such that the 
follower vehicles can coast. In this case, the coordination 
role of the platoon coordinator allows to avoid braking 
action to all vehicles. 

• Segment 2: Since no gear shift is simulated the vehicles 
are able to maintain the time gap requirement during the 
uphill. 

• Segment 3: due to the longer downhill compare to the 
first one, the platoon exhibits a different behavior. First, 
the platoon coordinator requires all vehicles to decrease 
the speed to the minimum allowed (in this simulation it 
is set to 19 m/s ) in order to hit the speed limit as late as 


Figure 16. Simulation results obtained using the proposed controller for a 
three-vehicle platoon while driving along the Sector A highlighted in Figure]!] 
The three vehicles are identical with parameters shown in Table]!] The first 
plot shows the road topography. For the explanation of the other plots refer 
to the caption of Figure |14| 

possible. Second, since the speed limit is reached despite 
the decrease of speed at the beginning of the downhill, the 
platoon coordinator requires the first vehicle to coast and 
the follower vehicles to brake slightly to maximize the 
efficiency. In fact, in this case, to require the first vehicle 
to fuel slightly and brake at the end of the downhill would 
be contradictory. 

In conclusion the platoon, under the control of the integrated 
control architecture, shows a fuel-efficient and smooth behav¬ 
ior. 

X. Conclusions and future works 
A. Conclusions 

In this paper we have presented a novel control architecture 
based on look-ahead control for fuel-efficient and safe HDV 
platooning. 

The use of a look-ahead control framework for HDV 
platooning has been first motivated by the analysis of real 
experiments. In particular in this analysis we concluded that 
the use topography information in order to predict the behavior 
of the vehicles and coordination between the vehicles can be 
beneficial for both fuel-efficiency and safety reasons. 

This led to the design of a novel control architecture for 
platooning. Such architecture is divided into two layers. A 
centralized higher layer, denoted as platoon coordinator, is 
responsible for the coordination of the platoon by defining 
a speed profile that is feasible and fuel-efficient for the entire 
platoon by exploiting preview topography information. Such 
speed profile is communicated to each block of the decentral¬ 
ized lower layer, denoted as vehicle controller layer. Within 
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each vehicle controller a model predictive control routine 
tracks the reference speed profile and generates the real-time 
desired acceleration for the low-level vehicle controller. 

The performance of such control architecture has been 
evaluated through the analysis of numerical experiments. In 
details, the performance of the two layers has been studied 
both separately and in conjunction. 

B. Future works 

In the modeling of the vehicle powertrain we have assumed 
that the gear ratio can be chosen on a continuous interval on 
a unlimited span. However, this is not typically the case in 
commercial HDVs, where usually the transmission is handled 
by a gearbox that introduces fixed gear ratios and power losses 
during the gear shifts. Therefore in some future works we 
want to investigate how the presence of a gearbox should 
be managed in a optimal way. The optimal engine speed 
as function of the generated power shown in Figure [6] and 
the knowledge of the current speed can be used to compute 
the instantaneous optimal gear ratio. However the power loss 
and the delay during the gear shift make the problem of 
when the HDVs should change gear (e.g., independently or 
simultaneously) and which gear they should engage not trivial. 

Secondly, we would like to investigate how external distur¬ 
bances, as traffic ahead or a vehicle cutting in the platoon, 
can be handled in an autonomous way within the platoon 
controller framework. So far, in fact, such disturbances have 
been assumed to be handled manually by the drivers. However 
the prediction of local traffic would allow the platoon to move 
fuel-efficiently and safely in it. 

Acknowledgment 

The authors gratefully acknowledge the European 
Union’s Seventh Framework Programme within the project 
COMPANION, the Swedish Research Council and the Knut 
and Alice Wallenberg Foundation for their financial support. 
The authors also thank Dr. Assad Alam and Dr. Henrik 
Pettersson from Scania for the fruitful discussions. 

References 

[1] Companion project, http://www.companion-project.eu/ Accessed: 2015- 
03-15. 

[2] Scania AB. Annual report 2013. 2013. 

[3] European Environment Agency. Specific C0 2 emissions per tonne-km 
and per mode of transport in Europe, 1995-2011. http://www.eea.europa. 
eu/data-and-maps/figures/ds_resolveuid/82SB924WG2 2012. Ac¬ 
cessed: 2015-04-16. 

[4] United States Environmental Protection Agency. EPA and NHTSA adopt 
first-ever program to reduce greenhouse gas emissions and improve fuel 
efficiency of medium- and heavy-duty vehicles. 2011. 

[5] A. Alam. Fuel-efficient Heavy-Duty Vehicle Platooning. Phd thesis, 
KTH Royal Institute of Technology, Stockholm, Sweden, 2014. 

[6] A. Alam, A. Gattami, and K.H. Johansson. An experimental study on the 
fuel reduction potential of heavy duty vehicle platooning. In Proceedings 
of the 13th International IEEE Conference on Intelligent Transportation 
Systems, Madeira, Portugal , pages 306-311, 2010. 

[7] A. Alam, J. Martensson, and K. H. Johansson. Experimental evalua¬ 
tion of decentralized cooperative cruise control for heavy-duty vehicle 
platooning. Control Engineering Practice, 38:11-25, May 2015. 

[8] A. Alam, J. Martensson, and K.H. Johansson. Look-ahead cruise 
control for heavy duty vehicle platooning. In Proceedings of the 16th 
International IEEE Conference on Intelligent Transportation Systems, 
The Hague, The Netherlands, pages 928-935, 2013. 


[9] R. Bellman. Dynamic programming. Princeton University Press, New 
Jersey, USA, 1957. 

[10] C. Bergenhem, Q. Huang, A. Benmimoun, and T. Robinson. Challenges 
of platooning on public motorways. In 17th World Congress on 
Intelligent Transport Systems, pages 1-12, 2010. 

[11] C. Bergenhem, S. Shladover, E. Coelingh, C. Englund, and S. Tsugawa. 
Overview of platooning systems. In Proceedings of the 19th ITS World 
Congress, Vienna, Austria, 2012. 

[12] F. Blanchini. Set invariance in control. Automatica, 35( 11): 1747—1767, 
1999. 

[13] C. Bonnet and H. Fritz. Fuel consumption reduction in a platoon: 
Experimental results with two electronically coupled trucks at close 
spacing. In Proceedings of the Future Transportation Technology 
Conference, Costa Mesa, USA, SAE Technical Paper 2000-01-3056, 
2010. 

[14] F. Borrelli, A. Bemporad, and M. Morari. Predictive control for 
linear and hybrid systems. In preparation, available online at 
http://www.mpc.berkeley.edu/mpc-course-material 2014. 

[15] F. Browand, J. McArthur, and C. Radovich. Fuel Saving Achieved in 
the Field Test of Two Tandem Trucks Fuel. Technical Report June, 
University of Southern California, 2004. 

[16] European Commission. White paper 2011 - Roadmap to a single 
European transport area: towards a competitive and resource efficient 
transport system. 2011. 

[17] European Commission. EU transport in figures - statistical pocketbook. 
2014. 

[18] W.B. Dunbar and R.M. Murray. Distributed receding horizon control for 
multi-vehicle formation stabilization. Automatica, 42(4):549-558, 2006. 

[19] International Transportation Forum. Transport outlook 2015. 2015. 

[20] E Hellstrom, J. Aslund, and L. Nielsen. Design of an efficient algorithm 
for fuel-optimal look-ahead control. Control Engineering Practice, 
18(11): 1318—1327, 2010. 

[21] E. Hellstrom, M. Ivarsson, J. Aslund, and L. Nielsen. Look-ahead control 
for heavy trucks to minimize trip time and fuel consumption. Control 
Engineering Practice, 17(2):245-254, 2009. 

[22] R. Horowitz and P. Varaiya. Control design of an automated highway 
system. Proceedings of the IEEE, 88(7):913-925, July 2000. 

[23] W.-H. Hucho, editor. Aerodynamics of road vehicles. Society of 
Automotive Engineers, USA, 4th edition, 1998. 

[24] K. H. Johansson, M. Tomgren, and L. Nielsen. Vehicle applications 
of controller area network. In Handbook of networked and embedded 
control systems, pages 741-765. Springer, 2005. 

[25] M. P. Lammert, A. Duran, J. Diez, K. Burton, and A. Nicholson. Effect 
of Platooning on Fuel Consumption of Class 8 Vehicles Over a Range 
of Speeds, Following Distances, and Mass. SAE International Journal 
of Commercial Vehicles, (30):7-9, September 2014. 

[26] J. Larson, K. Liang, and K. H. Johansson. A Distributed Framework 
for Coordinated Heavy-Duty Vehicle Platooning. IEEE Transactions on 
Intelligent Transportation Systems, 16(1):419—429, February 2015. 

[27] K. Liang. Coordination and Routing for Fuel-efficient Heavy-Duty 
Vehicle Platoon Formation. Licentiate thesis, KTH Royal Institute of 
Technology, Stockholm, Sweden, 2014. 

[28] G.J.L. Naus, R.P.A. Vugts, J. Ploeg, M.J.G. van de Molengraft, and 
M. Steinbuch. String-stable CACC design and experimental validation: 
a frequency-domain approach. IEEE Transactions on Vehicular Tech¬ 
nology, 59(9):4268—4279, 2010. 

[29] D. Norrby. A CFD study of the aerodynamic effects of platooning trucks. 
Master’s thesis, KTH Royal Institute of Technology, Stockholm, 2014. 

[30] L.E. Peppard. String stability of relative-motion PID vehicle control 
systems. IEEE Transactions on Automatic Control, 19(5):579—581, 
1974. 

[31] T. Sandberg. Heavy truck modeling for fuel consumption: Simulations 
and measurements. Technical report, 2001. LiU-TEK-LIC-2001:61, 
Thesis No. 924. 

[32] S.E. Shladover. PATH at 20. History and Major Milestones. IEEE Trans¬ 
actions on Intelligent Transportation Systems, 8(4):584-592, December 
2007. 

[33] S. Tsugawa. An overview on an automated truck platoon within the 
energy ITS project. In Advances in Automotive Control, volume 7, pages 
41^16, 2013. 

[34] V. Turri, B. Besselink, J. Martensson, and K. H. Johansson. Fuel- 
efficient heavy-duty vehicle platooning by look-ahead control. 53rd 
IEEE Conference on Decision and Control, pages 654-660, December 
2014. 


